MIT-CTP-2169 / ORNL-CCIP-92-15 

Kaon-Nucleon Scattering Amplitudes and Z*-Enhancements 
from Quark Born Diagrams 

T.Barnes 

Physics Division and Center for Computationally Intensive Physics 
Oak Ridge National Laboratory, Oak Ridge, TN 37831-6373 

and 

Department of Physics 
University of Tennessee, Knoxville, TN 37996-1200 

E.S.Swanson 
Center for Theoretical Physics 
Laboratory of Nuclear Science and Department of Physics 
Massachusetts Institute of Technology, Cambridge, MA 02139 



We derive closed form kaon-nucleon scattering amplitudes using the "quark Born 
diagram" formalism, which describes the scattering as a single interaction (here the 
OGE spin-spin term) followed by quark line rearrangement. The low energy 1=0 
and 1=1 S-wave KN phase shifts are in reasonably good agreement with experi- 
ment given conventional quark model parameters. For ki a b > 0.7 Gev however the 
1=1 elastic phase shift is larger than predicted by Gaussian wavefunctions, and we 
suggest possible reasons for this discrepancy. Equivalent low energy KN potentials 
for S-wave scattering are also derived. Finally we consider OGE forces in the re- 
lated channels KA, K*N and K*A, and determine which have attractive interactions 
and might therefore exhibit strong threshold enhancements or "Z*-molecule" meson- 
baryon bound states. We find that the minimum-spin, minimum-isospin channels 
and two additional K*A channels are most conducive to the formation of bound 
states. Related interesting topics for future experimental and theoretical studies of 
KN interactions are also discussed. 
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I. INTRODUCTION 



Kaon-nucleon collisions allow one to address many interesting problems in nuclear and 
hadron physics [|IJ. (By "kaons" we refer to the K + = us and K° = ds, generically K, as 
distinct from the K antikaons K~ = su and K° = sd.) Three familiar examples which we 
shall discuss below are 1) the origins of nonresonant "nuclear" forces in a system distinct from 
NN, 2) nuclear structure physics, using kaons as weakly scattered probes, and 3) searches 
for possible exotic Z* baryon resonances which couple directly to KN. More recently it has 
become clear that an understanding of KN scattering in nuclear matter is important in 
other areas, such as the interpretation of strangeness production in nuclear collisions and in 
two-kaon correlation measurements 0. 

Elastic KN scattering is a natural system for the study of nonresonant nuclear forces. 
Since the valence kaon wavefunction contains an s antiquark which cannot annihilate against 
the nonstrange nucleon state, direct production of conventional baryon resonances is ex- 
cluded. KN scattering is further simplified by the absence of one pion exchange, so one can 
study the nonresonant, non-OPE part of hadron scattering in relative isolation. Theoretical 
studies of KN nuclear forces are especially appropriate because there is already considerable 
experimental information on the elastic amplitudes and two-body inelastic reactions such as 
KN— > K*N and KN— > KA [T|,^-|^|. These experimental amplitudes provide stringent tests 
for models of hadronic interactions. The dominant S-wave elastic phase shifts are moder- 
ately well established, and the higher partial waves up to L=4 have been determined or 
estimated ||. The basic features of the elastic reaction are a strong repulsion in the 1=1 
S-wave, a weaker repulsion in the 1=0 S-wave, and an important spin-orbit interaction which 
is evident in the P-waves. The important low energy behavior of the 1=0 S-wave, in par- 
ticular the scattering length, is unfortunately not yet very well known. The experimental 
situation should improve considerably with the development of new hadronic facilities such 
as DA$NE and KAON 

KN scattering also has applications in nuclear physics; since the kaon-nucleon cross sec- 
tion is relatively small, kaon beams can be used as probes of nuclear structure. It would 
obviously be useful to understand the mechanism and properties of the kaon-nucleon in- 
teraction for this application. In view of this application one topic in this paper will be 
the derivation of effective low energy KN potentials from the nonrelativistic quark potential 
model. 

Another reason for interest in KN collisions is the possibility of producing flavor-exotic Z* 
baryon resonances. If discovered, these might be resonances with the quark valence structure 
q 4 s ||, where q = u or d. Such multiquark hadrons were widely predicted in the early days of 
the quark model ||, but it now appears that multiquark basis states usually do not support 
resonances, due to the "fall apart" effect |T0|JTT[] . The known exceptions are deuteronlike 
"molecule" states of hadron pairs, which should perhaps be classified as unusual nuclear 
species. (Nuclei themselves are excellent examples of the tendency of multiquark systems 
to separate into hadronic molecules.) In the meson- meson sector two KK molecule states 
are reasonably well established ||12|| , and there are several other meson-meson candidates 
T3fl . In the antikaon-nucleon sector the A(1405) is an obvious candidate KN molecule, and 



there presumably are other molecule states in channels with attractive interactions. Both 
the elastic reaction KN— >KN and inelastic processes such as KN— > K*N and KN— > KA can 
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be studied for evidence of exotic Z* baryon resonances. With a realistic model of hadronic 
interactions we might reasonably expect to predict the quantum numbers of exotic meson- 
baryon molecular bound states, should these exist. 

In this paper we apply the "quark Born diagram" formalism to KN scattering. In this 
approach we assume conventional nonrelativistic quark model wavefunctions for the asymp- 
totic hadrons, and calculate the Hamiltonian matrix element for scattering due to a single 
interaction between constituents in different incident hadrons. To form color singlet final 
states at lowest order one must then exchange constituents. The full Born amplitude is 
obtained by summing over all such processes coherently. (Similar constituent exchange 
mechanisms have been proposed for high energy hadron scattering [|14j], and there is strong 
experimental evidence in favor of this mechanism from large- 1 exclusive reactions This 
nonrelativistic Hamiltonian matrix element is then combined with relativistic phase space 
and kinematics to give results for differential cross sections, partial wave amplitudes and 
other scattering observables. In previous work we derived the elastic scattering amplitudes 
for 1=2 mi |16 |, 1=3/2 Ktt |17| and 1=1 KK |16[. (These cases were chosen because they are 



free of valence qq annihilation processes, which are known to be important if allowed.) We 
found good agreement with experimental titi and Ktt S-wave phase shifts given conventional 
quark model parameters. We have also applied similar techniques to pseudoscalar-vector 
and vector- vector meson channels |18| , and the results may have important implications for 
meson spectroscopy In Appendix C of [16|] we presented a diagrammatic representation 
of these techniques, with associated "Feynman rules" for the scattering diagrams. KN elastic 
scattering is also annihilation free and affords a nontrivial test of the quark Born formalism. 
KN elastic scattering has previously been the subject of numerous theoretical investi- 



gations. Meson exchange models have been applied in several studies 0], but these are 
difficult to justify fundamentally because the range of heavier meson exchange forces (~ 0.2 
fm) is much smaller than the minimum possible interhadron distance for two distinct hadrons 
(^ 1 fm) [ITJ. These models typically have many free parameters, which are not well estab- 
lished experimentally and are fitted to the data. Thus one is in effect simply parametrizing 
experiment. This type of model may be of theoretical interest as a parametrization of more 
fundamental scattering mechanisms which operate at the quark and gluon level, as it may 
be possible to relate the predictions of these different approaches. 

A quark and gluon approach to scattering using the P-matrix and bag model wave- 
functions was proposed by Jaffe and Low [pOj] . They suggested interpreting the multiquark 



clusters of the bag model not as resonances, but instead as the short distance parts of 
hadron-hadron scattering states. In principal this approach can be used to predict phase 
shifts, but in practice it has mainly been used to interpret experimental phase shifts in 
terms of P-matrix poles. This approach has been followed for KN by Roiesnal plf , who 
concluded that the KN data could indeed be interpreted in terms of poles approximately 
at the energies predicted by the bag model, but that the pole residues (coupling strengths 
to asymptotic KN channels) did not agree well with predictions. A more recent bag model 
calculation of KN scattering by Veit, Thomas and Jennings [^] used the cloudy bag model, 
which combines quark fields (in the baryon) with fundamental pseudoscalar meson fields 
in an effective lagrangian. This composite model leads to an 1=1 S-wave phase shift and 
a scattering length which are very similar to our result, but their 1=0 phase shift is much 
smaller than experiment. Although this cloudy bag approach gives promising numerical 
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results, it does not provide us with an understanding of the scattering mechanism at the 
quark and gluon level. 

Studies of the dominant S-wave KN scattering amplitudes in terms of quark model wave- 
functions and quark-gluon interactions have been published by Bender and Dosch [[^| (adi- 
abatic approach), Bender, Dosch, Pirner and Kruse (variational generator coordinate 
method, GCM) and Campbell and Robson (resonating group method, RGM). The 
large spin-orbit forces evident in the KN P-wave data have also been studied using similar 



quark model techniques, first qualitatively by Pirner and Povh |26| and later in detail by 
Mukhopadhyay and Pirner |27| (using GCM). The assumptions regarding dynamics, the 
scattering mechanism, quark model wavefunctions and the parameters used in these calcu- 
lations are very similar to our assumptions in this paper. The most important differences 
are that 1) our techniques are perturbative and allow analytic solution, and 2) we disagree 
about the size of the OGE contribution to KN scattering. Specifically, we find that OGE 
alone suffices to explain the observed 1=1 KN scattering length, whereas Bender et al. |2"4]| 
conclude that OGE is too small, and that a Pauli blocking effect is dominant in 1=1. Camp- 
bell and Robson |25[ similarly found that the experimental 1=1 phase shift was larger than 
their predictions, which were based on generalizations of Gaussian wavefunctions and a full 
OGE and confining interaction. 



II. CALCULATION OF KN AND RELATED SCATTERING AMPLITUDES 



a) Hamiltonian and hadron states 

Our technique involves a Born order calculation of the matrix element of the Hamiltonian 
between asymptotic hadron states in the nonrelativistic quark model. In the KN case the 
dominant interaction was previously found by Bender et al. p|| to be the spin-spin "color 
hyperfine" term. A similar conclusion has been reached for the NN interaction [TT|,|2S|. Here 
we shall adopt this approximation and neglect the other OGE and confining terms. Thus, 
our scattering amplitude is proportional to the matrix element of 
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between asymptotic KN states. {Tf is the color matrix for quark or antiquark i, which is 
A a /2 for quarks and — (X T ) a /2 for antiquarks.) Although we shall quote results for arbi- 
trary asymptotic hadron wavefunctions, we shall specialize to Gaussian wavefunctions for 
our numerical results, as these allow closed form derivation of scattering observables. Our 
momentum space Gaussian wavefunctions for the kaon and nucleon are conventional quark 
model forms, 



Ykaon [Prel] 
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and 



0nucleon(pi,P2,P3) = eX P j — j . (4) 

The parameters a and /? are typically found to be ~ 0.3 — 0.4 Gev in hadron phenomenology. 
These are relative momentum wavef unctions, and have an implicit constraint that the con- 
stituent momenta add to the hadron momentum. In the full momentum space wavefunction 
there is an overall delta function that imposes this constraint; 

®kaon{Pq,Pq] Ptot) = 4>kaon{Prel) S(P tot ~Pq~ Pq) , (5) 
®nucleon(Pl,P2,P3] Ptot) = <t>nucleon{Pl, P2, Pz) S{P to t ~ Pi ~ P2 ~ Pi) ■ (6) 

The normalizations are 

( $ kaon ( Plot ) I &kaon ( Ptot ) ) 



and 



dpdpdp' dp' % aon (p\ p'; PU)$kaon(p, P; Ptot) = $(Ptot ~ PL) (7) 

*&nucleon(Ptot) l^nucleon (Pot)) 

dpi dp 2 dp 3 dpi ' p 2 ' dp 3 ' KucleoniPl P2 ', P3 '] Ptot ) ^nucleon (Pi j P2,P3]Ptot) 



= S(P tot - PL) . (8) 
Since these state normalizations are identical to those used in our previous study of K7r 



scattering ]I7) we can use the relations between amplitudes and scattering observables given 
in that reference. 

The color wavefunctions for the asymptotic hadrons are the familiar color singlet states 

\meson) = —= 6„ \vi) (9) 



1,1- 



^1,3 



and 



\baryon) = ^2 A= e ijk \ijk) . (10) 

i,j,k=l,3 V 

Our spin-flavor states for the meson and baryon are standard SU(6) states, but we have 
found it convenient to write the baryon states in an unconventional manner, as the usual 
quark model conventions are unwieldy for our purposes. First, to establish our notation, the 
spin-flavor K + state is 
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|K+) = i(|u + 5_)-|u_5 + )) • (11) 

For quark model baryon states it is conventional to assign each quark a fixed location in the 
state vector, as though identical quarks were distinguishable fermions. One then explicitly 
symmetrizes this state. Thus for example one writes the normalized A + (S Z = +3/2) state 
as 

| A+(+3/2)> = -y= (\u + u+d+) + \u + d + u + ) + \d + u + u + )^ (12) 

and treats each basis state as orthogonal. Note however that this is not the usual way to 
represent multifermion states. In standard field theoretic usage each of these basis states is 
identical to the others, to within an overall phase. In this language the normalized A + (+3/2) 
state is simply 

|A+(+3/2)> = -J= \u + u + d + ) , (13) 

which we could equally well write as \u + d + u + ) / \[2 or \d + u + u + ) / v2. The advantage of using 
field theory conventions becomes clear in calculating nucleon matrix elements. For example, 
the usual quark model proton state is 



|P(+l/2))= (^2|w_ 




h cL>- 


- \u + u. 


-d+) ~ 1 


+ 2|m_ 


h cL 


-«+) " 


- \u + d A 


-«-> " 1 


+ 2\d. 




-«+) - 


- \d + u^ 


M-) — 



(14) 



and in comparison this state in field theory conventions is 




|P(+l/2)> = f- {^f-} - f 3 \u + u. d+ ) . (15) 

Use of the latter form, with all permutations of quark entries allowed in matrix elements, 
reduces the number of P^P terms from 81 (many of which are zero) to 4. Of course the 
results are identical, as these are just different conventions for the same state. 

b ) Enumeration of quark line diagrams for KN 



Now we consider KN scattering. As explained in reference [fLql , we begin by determining 
the matrix element of the scattering Hamiltonian (1). First we factor out the overall mo- 
mentum conserving delta function and then derive the remaining matrix element, which we 
call hfi] 

f(KN\H scat \KN}i = h fi 6(P f - Pi) . (16) 

We will discuss one part of the calculation in detail to explain the techniques, and then 
simply quote the full result. Specializing to the spin up 1=1 case K + P(+l/2)^K + P(+l/2), 
we require the matrix element of the scattering hamiltonian (1) between initial and final K + P 
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states with color and spin-flavor wavefunctions given by (9,10) and (11,15) respectively. Since 
the kaon and proton states (11) and (15) are each the sum of two terms, the full amplitude 
for K + P(+l/2)— >K + P(+l/2) is a weighted sum of 16 subamplitudes. We shall consider the 
subamplitude for \u + sJ){\u + u + d_)/\/2} — > \u + sJ)\u + U-d+) , which we call hff', in detail 
for illustration. 

We begin by constructing all allowed quark line diagrams and their associated combi- 
natoric factors. First we arrange the initial and final states with their normalizations on a 
generic scattering diagram, 



u + 

(17) 

u + 
d + 



Now we connect the initial and final lines in all possible ways consistent with flavor 
conservation. For the d quark and the s antiquark this choice is unique. For the final 
meson's u quark however there are two choices for which initial baryon's quark it originates 
from. Similarly the initial meson's u quark can attach to either of two final baryon u quarks. 
Thus we have four quark line diagrams. We may immediately simplify the diagrams; since 
the baryon wavefunctions are symmetric, we may permute any two initial or final baryon 
lines and obtain an equivalent diagram. We use this symmetry to reduce all diagrams to a 
"standard form" in which only the meson's quark and the upper baryon quark are exchanged. 
The two choices for the initial baryon's spin up u + quark are thus equivalent, and contribute 
an overall combinatoric factor of two. The final baryon's quarks however give inequivalent 
diagrams, one being nonflip (u + (K) — > u + (P)) and the other spin flip (u + (K) — > u_(P)). 
(No polarization selection rules are being imposed yet, only flavor conservation.) Thus our 
amplitude leads to the line diagrams 



h e.g. 



1 

V2 



U + 
S- 



U+ 

u + 

d^ 



We next "decorate" each of these line diagrams with all possible single interactions (1) 
between a quark (or antiquark) in the initial meson and a quark in the initial baryon. There 
are six of these per line diagram (two choices in the meson times three in the baryon), so we 
have a total of twelve scattering diagrams to evaluate. However in this case all but one are 
trivially zero. Note that in the first line diagram we must flip the spins of u and d quarks in 
the initial baryon to have a nonzero contribution. This however is not part of our scattering 
interaction, which operates between pairs of constituents in different initial hadrons. The 
Si ■ Sj interaction either flips a pair of spins in different incident hadrons (through S + S- 
or SS+ terms) or leaves all spins unchanged (through S Z S Z ). Thus, the transition in the 

— * — * 

first line diagram cannot occur through a single Si ■ Sj interaction. For the second diagram 
however there is a single nonvanishing transition, in which the initial meson's u + quark and 
the baryon's gL quark interact through the spin flip operator; 




c) Independent quark and gluon diagrams and their spin and color factors 

Finally we require the spin, color, overall phase and spatial factors associated with this 
and the other independent diagrams. There are only four independent quark and gluon 
diagrams, since all others can be obtained from these by permutation of lines. These four 
diagrams are 
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The spin factor is simply the matrix element of Si ■ Sj for scattering constituents i and j, 
evaluated between the initial and final (qq)(qqq) spin states. This is 1/2 if both spins i and 
j are antialigned and both flip, +1/4 if the spins are aligned and neither flips, and —1/4 if 
they are antialigned and neither flips. All other cases give zero. All spectator spins must 
not flip or the overall spin factor is trivially zero. 

The color factor can be evaluated using the states (9), (10) and standard trace techniques, 



as in (51) of reference [16]. The result for each diagram is 

/color (£>i) = +4/9 , (24) 



^color (D 2 ) = 


-2/9 , 


(25) 


/color (-D3) — 


-4/9, 


(26) 


/color (-D4) = 


+2/9 . 


(27) 
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d) "Diagram weights" for KN scattering 

We conventionally write the meson-baryon hfi matrix elements as row vectors which 
display the numerical coefficient of each diagram's spatial overlap integral. Thus, 



h 



fi 



Wl, U>2, W 3 , W 4 



represents 



hfi — W\ -^space(-^l) "I" ^2 -^space(-02) "I" ^3 -^spacc(-Os) ^4 -^space 



(28) 



(29) 



This notation is useful because the diagram weights {wi} are group theoretic numbers that 
obey certain symmetries, whereas the spatial overlap integrals are complicated functions that 
depend on the specific spatial wavefunctions rather than the symmetries of the problem. 



As an illustration, our practice subamplitude h e jf' is 



h% 8 ' = V2 •(£)■(- |)-WA) 



2 / V 9, 

(using the spin and color matrix elements given above), which we abbreviate as 



ue.g. 



o, - 



x/2 
9 ' 



0, 



(30) 



(31) 



This completes our detailed derivation of h e jf for the subprocess \u+S-){\u + u + d_) / f y2} — > 
\u + s_)\u + u_d + ). 

Proceeding similarly, we have derived the weights for the full KN elastic scattering am- 
plitudes, given the states (11) and (15) and their isospin partners. These are 



h™(I = 0) 



and 



1=1) 



0, -, 0, - 

' 6' ' 6 



1111 

3' 18' 3' 18 



(32) 



(33) 



For numerical estimates of these amplitudes we require the spatial overlap integrals, which 
we shall evaluate explicitly with Gaussian wavefunctions. 

e) Spatial overlap integrals 

The spatial overlap integrals represented by the four diagrams D\ . . . D4 may be deter- 
mined using the diagrammatic techniques discussed in Appendix C of reference ||161 . These 
are formally 30-dimensional overlap integrals (three dimensions times ten external lines), 
but twelve integrations are eliminated by external momentum constraints and an additional 
nine are eliminated by the unscattered spectator lines. This leaves a nontrivial 9-dimensional 
overlap integral for each diagram. We give the initial meson a label A, with three-momentum 
also called A and quark three-momentum a and antiquark momentum a, and we similarly 
label the initial baryon B, the final meson C and the final baryon D. Since we choose to 
evaluate these integrals in the CM frame we use the momentum substitutions B = —A and 
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D = —C. We also introduce a nonstrange to strange quark mass ratio p = m q /m s . With 
these substitutions the four spatial overlap integrals are 

I sP ace{Dx) = HI dddb, db 2 <\> A {2a - <f>* c (2a + ^- - 2A) 



(h, b 2 , -A-h- h) <f>* D {bi +A-C, b 2 , -A-h- h) , (34) 



Is P aee(D 2 ) = +^ — ^ 1 1 1 db, dcdd, <\> A (2 C - - 2C) 0* C (2 C ^ ' 



3m2 {2-KfJJJ 1 1 1 + p ^ cv l + p J 

■ (j) B (h, c, -A — bi — c) <ff D (di, A-C + b 1 + c-d 1 ,-A-bi-c) , (35) 

8-Ka s 1 f r f lT * , , 2pA 2pC . 

WAO = +^PJ^JJJ dadb 2 dc M2a - —) 4>c(2c-—) 



{a-A + C, b 2 , -a-h- C) <j>* D {a, b 2 , -a - b 2 - C) , (36) 



8fca s 1 f f f 7 _ ,r* 2pA 2pC 

I sp ace{D 4 ) = + ^-^-p yr / / / aaci&i dc 0^(20 - — — ) c (2c 



3m^(2vr) 3 yyy 1 ^ v l + p /Y ^ v l + p y 

• 0b(6i, c, -A - 61 - c) 0o(A - C - a + h + c, a, -A - h - c) . (37) 

There are many equivalent ways to write these integrals which arise from different choices 
of the variables eliminated by momentum constraints. 

Note that the overall coefficients of these integrals are positive, although the coefficient 
of H scat (1) is negative. This is because there is an overall phase factor of (—1) for each 
diagram Di . . . D4, due to anticommutation of quark creation and annihilation operators in 
the matrix element. Here we incorporate this phase, which we call the "signature" of the 
diagram [16], in the spatial overlap integrals. The signature is equal to (—l) Nx , where N x is 
the number of fermion line crossings. For diagrams D% . . . D4 above N x = 3, so the signature 
is 



' signature 



(38) 



Note that a diagram in nonstandard form, such as the kaon's quark line crossing to the 
second baryon quark, can have a (+1) signature; in the full hfi matrix element this is 
compensated by a change in sign of the color factor. 

We explicitly evaluate these overlap integrals using the Gaussian wavefunctions (2) and 
(4). For Gaussians the integrals factor into products of three 3-dimensional integrals, and 
the results are all of the form 



8rra s 1 f 2 1 
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where P 2 m is the modulus of each hadron's three-momentum in the CM frame, p = cos(6cm) 
where $cm is the CM scattering angle, and the constants rji, and Bi are functions of a, (3 
and p. Bi > implies forward peaked scattering and Bi < is backward peaking. The pure 
exponential dependence in P 2 m and P 2 m p is a consequence of the Gaussian wavefunctions 
and the contact interaction. Introducing the ratio g = (a/j3) 2 , these constants are 

Vi = 1 (40) 



_ 2p 2 + 4p+(30 + 2) 
1 6(l + p) 2 a 2 1 ' 



B 1 = A 1 , (42) 



^ = I TiT-Q) (43) 



_ (400 + 3)p 2 + (320 + 6)p + (2l0 2 + 28g + 3) 
2 6(70 + 6)(l + p) 2 a 2 1 ' 



B (-8g + l)p 2 + 2p+(7g 2 + 8g + l) 
2 2{7g + 6){l + p) 2 a 2 ' 1 } 



6 \ 3 / 2 

* = P( 5T 3) (^) 



A (10g + Q)p 2 + (8g + 12)p+(7g + 6) 
3 6((? + 3)(l+p) 2 a 2 1 ] 



_ (-0 + l)p 2 + 2p+(0 + l) 
* S " (0 + 3)(l + p) 2 a 2 ' (48) 



/ 120 \ 3 / 2 

r?4 - p l(20 + 3)(0 + 2)J 



(49) 



(2O0 2 + 400 + 3)p 2 + (40 2 + 140 + 6)p + (5g 2 + lOg + 3) 

6(20 + 3)(0 + 2)(l + p) 2 a 2 ' [M) 



_ (-40 2 - 80 + l)p 2 + (-40 2 - 60 + 2)p + (0 2 + 20 + 1) 

2(20 + 3)(0 + 2)(l + p) 2 a 2 • 1 ' 
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These results were derived at MIT and ORNL independently using MAPLE and MAC- 
SYMA algebra programs respectively. 

Some important properties of these diagrams, specifically which are forward peaked or 
backward peaked processes, and which diagrams dominate at high energies, can be inferred 
by inspection. The leading diagram in the high energy limit is D\, which is a forward peaked 
exponential in t. The other diagrams are exponentially suppressed in s and are also forward 
peaked, with the single exception of D4. Note that for plausible values of g ~ 1 and p ~ 0.6 
this diagram leads to a backwards peak (i? 4 < 0). These properties have a simple common 
origin; since we are scattering through a hard delta function interaction, the only angular 
dependence comes from overlap suppression due to the spectator lines. A spectator line 
which is required to cross into the other hadron gives an especially large suppression at 
high energies and small angles. The amplitude for a crossing spectator line is maximum for 
backscattered hadrons; in this case the crossing spectator is actually continuing to move in 
a new hadron with the same momentum vector as the hadron it originally resided in. 

The first diagram Di has no crossing spectators, so it is not suppressed in s; only the hard 
scattered constituents are required to cross into different hadrons. In diagrams D2 and D3 
one spectator line is required to cross to a different hadron, so there is some suppression with 
increasing s. Since two spectators do not cross, they dominate the angular dependence, and 
the scattering is forward peaked. Diagram D4, the backward peaking process, is qualitatively 
different because two spectator lines are required to change hadrons, and only one spectator 
does not cross. In this case "backwards" meson-baryon scattering actually corresponds to 
forward scattering for the two crossing spectator lines, which is obviously preferred. This 
description attributes backward peaks, which might otherwise appear counterintuitive, to 
the obvious mechanism of "minimum spectator suppression" at the quark level. 

f) KN phase shifts and scattering lengths 

Given the diagram weights (32-33) and our results (40-51) for the Gaussian overlap 
integrals, we have completed the derivation of the Hamiltonian matrix element hfi for KN 
elastic scattering. Since we have used the same normalization for KN states as in our previous 
discussion of Kir scattering [17] we can use the same relations derived there to relate hfi to 



scattering variables. First we consider the elastic phase shifts, which are given by 

= jj^^eEe f h ™ Plifl)dfl . (52) 

(E K + E N ) J-i 

Using the integral j\ e bu 'Pe(n)d/j, = 2^(6), we find 

5fN = _^ P^EkE n * ^ ^ eM _ AiP 2j m p L) (53) 
om, \&k + En) i=1 

where one specifies the isospin state 1=0 or 1=1 through the choice of the diagram weights 

As we approach the KN threshold the S-wave phase shift is asymptotically linear in P cm , 
and the coefficient is the scattering length aj. Since the exponential and the %q Bessel function 
are both unity in this limit, we recover a relatively simple result for the KN scattering length, 

KN 4a, M K M N ' 
<" = ~MAM K + ■ (54) 
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Since the coefficients {n} are relatively simple functions, we can write these scattering 
lengths as simple functions of a s /m 2 , p = m q /m s , the meson-baryon relative scale parameter 
g = (a//?) 2 and the physical masses M K and M N . The results are 

KN = 4a a M K M N 
ai=1 ~ 3w? q {M K + M N ) 



and 



1 1/125 
3 + 18 V7# + 6 



3/2 ]_ 
+ 3 P 



6 \ 3 / 2 J_ 
5 + 3J + 18 



12 5 



\ 3 / 2 



(2g + 3)(g + 2). 



(55) 



l 7=0 



4a s M K M, 



N 



3ml (M K + M N ) 



M / 12g \ 3 / 2 1 / 12g y/a 
.6V7^ + 6y 6 P V(2^ + 3)(^ + 2)y 

The basic features of the low energy KN interaction, a repulsive 1=1 S-wave and a repul- 
sive but less strong 1=0 S-wave, are already evident in these formulas. (The parameter 
g = (a/(3) 2 is constrained by quark model phenomenology to be comparable to unity.) De- 
tailed numerical results for the scattering lengths and phase shifts and a comparison with 
experiment are presented in the next section. 



(56) 



III. COMPARISON WITH EXPERIMENT 

a) Scattering lengths 

Before we discuss our numerical predictions we first review the status of the experimental 
scattering lengths. Since there are unresolved disagreements between analyses in the 1=0 
channel, we have compiled relatively recent (since 1980) single-energy S-wave phase shifts 
for our discussion. These are in chronological order Martin and Oades |3(J (Aarhus and UC 
London, 1980); Watts et al. (QMC and RAL, 1980); Hashimoto [|32J (Kyoto and VPI, 
1984); and Hyslop et al. § (VPI, 1992). The 1=1 data set analysed by Arndt and Roper |33 



(VPI, 1985) was incorporated in the 1992 VPI simultaneous analysis of 1=0 and 1=1 data, 
so we shall not consider it separately. The energy dependent parametrizations of Corden et 
al. |34j] and Nakajima et al. |3l| are not included in our discussion. 



In Fig.l we show these experimental 1=0 and 1=1 S-wave phase shifts versus P cm = |P cm |. 
The linear low energy behavior which determines the scattering length is evident in the 
1=1 data, and Hyslop et al. cite a fitted value of afJi = —0.33 fm. Previous analyses 
(summarized in [||] and H) have given values between —0.28(6) fm |36f and —0.33 fm |5|.|37|. 
A more useful way to present the S-wave phase shift data is to display 5q/ P cm versus P 2 m ; 
the intercept is the scattering length, and the slope at intercept determines the effective 
range. In Fig. 2 we show the S-wave phase shifts in this manner; an 1=1 scattering length 
of about —0.31(1) fm is indeed evident, which we shall take as our estimated experimental 
value. 
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Unfortunately the 1=0 scattering length is much less well determined, as is evident in 
Figs.l and 2. Previous (favored) solutions up to 1982 are summarized in Table 2.3 of 
n, and range between +0.02 fm and — O.lli^Ew There appear to be two sets of low 
energy values in the data of Fig.l, a smaller phase shift from the Aarhus-UCL and QMC- 
RAL collaborations and a larger one from from the Kyoto- VPI and VPI analyses. Below 
P cm = 0.4 Gev the Kyoto- VPI and VPI results are larger than Aarhus-UCL and QMC-RAL 
by about a factor of two. The VPI group actually cite a scattering length of afj^ = 0.0 
fm, although this requires rapid low energy variation below the first experimental point 
(compare their Fig.l(a) with the 1=1 phase shift in their Fig.2(a), which is constrained by 
experiment at lower energy and shows the expected ^/T[ ab oc P cm S-wave dependence). Since 
the 1=1 phase shift is close to linear for P cm < 0.4 Gev {ki a b < 0.7 Gev), we will assume 
that the zero 1=0 scattering length quoted in [E) is an artifact of their fit, and that the 
actual 1=0 phase shift is approximately linear in P cm for P cm < 0.4 Gev. We can then read 
the 1=0 scattering length from the intercept in Fig. 2. From the figure we see that a naive 
extrapolation to threshold leads to scattering lengths of about —0.09(1) fm and —0.17(2) 
fm respectively from the two sets of references. In summary, the experimental phase shifts 
shown in Fig. 2 suggest to us the scattering lengths 



\(expt. 



-0.31(1) fm 



afj^(expt.) 



Aarhus-QMC-RAL—UCL 



Kyoto-VP I 



-0.09(1) fm 



-0.17(2) fm . 



(57) 



We emphasize that the 1=0 values are our interpretation of the data from Fig. 2, and the 
references cited quote smaller scattering lengths that we believe the data does not support. 
As the values of the 1=0 scattering length and low energy phase shifts are controversial, an 
accurate determination should be a first priority at a kaon facility. 

To compare our predictions with experiment we first use a "reference parameter set" 
with conventional quark model parameters. The hyperfine strength is taken to be a s /m 2 q = 
0.6/(0.33) 2 Gev -2 , and the nonstrange to strange quark mass ratio is p = m q /m s = 
0.33 Gev /0.55 Gev = 0.6. The remaining parameter in the scattering length formulas 
is g = (a/f3) 2 , the ratio of baryon to meson width parameters squared. These parame- 
ters are rather less well determined phenomenologically. For baryons, values in the range 
a = 0.25 — 0.41 Gev have been used in nonrelativistic quark model studies [38-iTJ]. Isgur 
and Karl j|T] originally used a = 0.32 Gev for spectroscopy, but Copley, Karl and Obryk 



42[ had earlier found that the photocouplings of baryon resonances required a somewhat 
0.41 Gev, which may be a more realistic estimate 



larger value of a = 0.41 Gev, which may be a more realistic estimate P8"| , |3"9"|] because it is 
less sensitive to short distance hyperfine matrix elements. This larger value was also found 
by Koniuk and Isgur [43 for baryon electromagnetic transition amplitudes. Here we take 
a = 0.4 Gev as our reference value. For mesons, studies of various matrix elements have led 
to values of (3 = 0.2 — 0.4 Gev In our previous study of 1=2 nir scattering we found a 
best fit to the S-wave phase shift data with j3 = 0.337 Gev. Here we use a similar (3 = 0.35 
Gev as our reference value; if the quark Born formalism is realistic we should use essentially 
the same meson parameters in all reactions. 
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With our reference parameter set and physical masses Mk = 0.495 Gev and = 0.940 
Gev, our formulas (55) and (56) give 

af=i(ref. set) = -0.35 fm (58) 

and 

af3>(ref. set) = -0.12 fm . (59) 



In view of our approximations, the parameter uncertainties, and the uncertainties in 
the 1=0 data, these scattering lengths compare rather well with experiment. Note that our 



conclusions differ from those of Bender et al. [ 24"| , who reported that the OGE contribution 



to 1=1 scattering was too small to explain the observed S-wave phase shift. We discuss this 
disagreement further in the appendix. 

Now suppose we attempt to fit our estimated experimental values of the scattering lengths 
(57) by varying our parameter set. It is useful to fit the ratio afj^/afj^, since this involves 
only p and the width parameter g. We have fixed p = 0.6, and in any case we find that 
af^/af^ is insensitive to p, so only g remains as an important parameter. In Fig. 3 we show 
the predicted ratio of KN scattering lengths as a function of a/ (3. The two experimental 
ratios assuming the values in (57) are also indicated. The larger ratio af= Q /afI{ = 0.17/0.31 
requires a//3 = 1.91, rather far from typical quark model values. Fitting the smaller ratio 
a 7=o/ a 7=i = 0.09/0.31 requires a/(3 — 1.02, which is more representative of quark model 
parameters. An accurate determination of the 1=0 KN scattering length through direct low 
energy measurements, rather than by extrapolation, would be a very useful experimental 
contribution; this would allow a more confident test of our results and those of other models 
(as shown for example in Table 6-4 of Hyslop 



b) S-wave phase shifts 

The S-wave KN phase shifts predicted by (53) with I = given the "reference parameter 
set" a s /m 2 q = 0.6/(0.33) 2 Gev" 2 , p = m q /m s = 0.6, a = 0.4 Gev and (3 = 0.35 Gev are 
shown as dashed lines in Fig. 4. As we noted in the previous section, this parameter set 
gives reasonable scattering lengths, although the 1=0 scattering length is not yet very well 
established experimentally. 

At higher energies the reference parameter set predicts an 1=1 phase shift that retreats 
more quickly with energy than is observed experimentally; in Fig. 4 we see a rapid departure 
of theory and experiment above P cm = 0.4 Gev {ki^ = 0.7 Gev). This is near the opening 
of the inelastic channels KA and K*N, as indicated in Fig. 4. Two possible reasons for 
this discrepancy are 1) inelastic effects of the channels KA, K*N and K*A, which should 
become important just where theory and experiment part, and 2) short distance components 
in the meson and baryon wavefunctions that are underestimated by the smooth Gaussian 
wavefunctions (2) and (4). 

Although inelastic effects are certainly important experimentally |||4[], the most impor- 
tant low energy inelastic process is P-wave KA production Qj. Hyslop et al. |JJ similarly 
find relatively small inelasticities in the 1=1 KN S-wave, with 7] > 0.9 for P cm < 0.68 Gev. 
At the end of this range our predicted phase shift given the reference parameter set is only 
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about half the observed value, so it appears unlikely that the discrepancy is mainly due to 
inelastic channels. 

A second possible reason for the discrepancy is a departure of the hadron wavefunctions 
from the assumed single Gaussian forms at short distances; both the meson qq states and the 
baryon qq substates experience attractive short distance interactions from the color Coulomb 
and hyperfine terms (for spin singlets), which will lead to enhancements of the short distance 
components of their wavefunctions and increased high energy scattering amplitudes. If this 
is the principal reason for the discrepancy, we would expect a global fit to the S-wave phase 
shifts to prefer a smaller hadron length scale. In Fig. 4 we show the result of a three- 
parameter fit to the full 1992 VPI 1=0,1 energy independent S-wave data set |J, letting 
a s /m 2 ql a and (3 vary and holding p = 0.6 fixed. The fit is shown as solid lines, and is 
evidently quite reasonable both near threshold and at higher energies. The fitted parameters 
are a s /m 2 q = 0. 59/(0. 33) 2 Gev~ 2 , a = 0.68 Gev and (3 = 0.43 Gev; the hyperfme strength 
is a typical quark model result but the width parameters a and (3 are about 1.5 times 
the usual quark model values. Thus, a fit to the S-wave VPI data using single Gaussian 
wavefunctions requires a hadron length scale about 0.7 times the usual scale. This result is 
largely independent of the data set chosen, since it is driven by the large 1=1 phase shift, 
which shows little variation between analyses. Evidently the predicted S-wave phase shifts 
at higher energies are indeed very sensitive to the short distance parts of the wavefunction; 
this supports our conjecture that the discrepancy at higher energies is an artifact of our 
single Gaussian wavefunctions. A calculation of these S-wave phase shifts using realistic 
Coulomb plus linear plus hyperfine wavefunctions is planned |3£| and should provide a very 
interesting test of the quark Born formalism. 

c) Higher-L partial waves, spin-orbit and inelastic effects 

In addition to the S-wave phase shifts, higher-L KN elastic phase shifts and properties 
of the inelastic reactions KN— »K*N, KN— >KA and KN— >K*A have been the subjects of 
experimental investigations. These studies have found important effects in the L>0 partial 
waves which are beyond the scope of the present paper. 

One especially interesting effect is a remarkably large spin-orbit interaction in the 1=0 KN 
system; the L=l states have a large, positive phase shift for J=l/2 and a weaker, negative 
phase shift for J=3/2 (see || and references cited therein). This spin-orbit interaction 
cannot arise in our quark Born amplitudes given the approximations we have made in this 
paper; since we have incorporated only the spin-spin hyperfine interaction in single hadronic 
channels, our phase shifts (53) are functions of the total hadronic L and S but not J. Some 
but not all of this spin-orbit interaction may simply require incorporation of the OGE spin- 
orbit term; Mukhopadhyay and Pirner [27] found that the quark spin-orbit interaction was 
sufficient to explain the sign and magnitude of some of the weaker KN spin-orbit forces, 
but that the 1=0, J=l/2 phase shift was much too large to be explained as an OGE force. 
The strong KN spin-orbit forces might conceivably be due to couplings to inelastic channels; 
since the available mixing states and their couplings to KN are J-dependent, they might 
lead to effective spin-orbit forces at the hadronic level, even if we do not include spin-orbit 
forces at the quark level. We hope to treat this interesting possibility in a future study of 
coupled channel effects using the quark Born formalism. 

Since we do not have a model of the large spin-orbit effect it is not appropriate to include 
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a detailed discussion of our predicted amplitudes and cross sections at higher energies, where 
higher partial waves are important. In the interest of completeness however we will briefly 
discuss our predicted differential cross section at high energy, since we previously noted that 
we found an exponential in t in 1=2 tttt scattering, reminiscent of diffraction in magnitude 
but not in phase ]16[. The differential cross section in this unequal mass case is related to 
the hfi matrix element (29) by 



da - 

— = 47T 5 

dt 



s 2 - (M% - Ml 



2 \2 



(s - (M N + M K f)(s - (M N - M K f) 



\h 



fi\ 



(60) 



For KN scattering in the high energy limit only the contribution from diagram Di (20) 
survives, and we find 



(61) 



Thus we again find an exponential in t at high energy, with a slope parameter (41) that is 
numerically equal to 



b = A 1 = 3.7 Gev" 



(62) 



given our reference parameter set. This is similar to the observed diffractive 1=1 KN slope 



parameter |44| of 



b(expt.) ps 5.5 — 5.9 Gev" 



(63) 



The normalizations of the theoretical and experimental 1=1 high energy differential cross 
sections however differ by about an order of magnitude, and are 1.8 mb Gev™ 2 (reference 
parameter set) versus ps 15 mb Gev -2 (experiment [f|4]]). We noted a similar tendency for the 
reference parameter set to underestimate high energy amplitudes in our discussion of the S- 
wave phase shifts, which we attributed to the single Gaussian wavefunction approximation. 
One interesting prediction is that 1=0 KN scattering should have no diffractive peak in 
the high energy limit, since it has has W\ = 0; unfortunately there is no 1=0 high energy 
data to compare this prediction with. A serious comparison with high energy scattering will 
presumably require the use of wavefunctions with more realistic high momentum components 
as well as the incorporation of inelastic channels, which may strongly affect the elastic 
amplitudes. 



IV. KN EQUIVALENT POTENTIALS 

Sufficiently close to threshold our quark Born scattering amplitudes can be approxi- 
mated by local potentials. These potentials are useful in applications such as multichannel 
scattering and investigations of possible bound states, which are easiest to model using a 
Schrodinger equation formalism with local potentials. There are many ways to define an 
equivalent low energy potential from a scattering amplitude such as hj^ several such pro- 
cedures are discussed in HHH^] and in Appendix E of Ul6"f . Of course effective potentials 



18 



extracted using different definitions can appear to be very different functions of r although 
they lead to similar low energy scattering amplitudes. 

One approach to defining an equivalent potential is to derive a potential operator V op (r) 
which give the scattering amplitude hfi in Born approximation. This "Born-equivalent po- 
tential" technique is discussed in reference [|5] and in Appendix E of flBfl ; it has been tested 



on the OGE interaction, from which one recovers the correct Breit-Fermi Hamiltonian at 
0(v 2 /c 2 ) ||45||. To derive the Born-equivalent potential we reexpress our scattering amplitude 



in the CM frame as a function of the transferred three-momentum q = C — A and an orthog- 
onal variable V = (A + C)/2. We then expand the scattering amplitude in a power series in 
V and equate the expansion to the Born expression for nonrelativistic potential scattering 
through a general potential operator V op (r), which may contain gradient operators. The 
leading term, of order V°, gives the Born-equivalent local potential V(r). 

In this meson-baryon scattering problem our Hamiltonian matrix elements are of the 
form 



8na s 1 

h fi 



j^WiVi exp {-(A- B lf i)P 2 \ . (64) 
i=i L ' 



3m 2 (2tt) 3 

Making the required substitutions P 2 m = V 2 + <f 2 /4 and P 2 m p = V 2 — <f 2 /4 and Fourier 



transforming with respect to q as in |16| gives the equivalent low energy KN potential 

VMr) - ^ g JX ^ m exp { - rV (Ai + B,)}. (65) 

Thus our Born-equivalent meson-baryon potentials are sums of four Gaussians, one from 
each inequivalent quark Born diagram, weighted by the diagram weights of that channel. 

The potentials for 1=0 and 1=1 with our reference parameter set a s = 0.6, m q = 0.33 
Gev, p = m q /m s = 0.6, a = 0.40 Gev and [3 = 0.35 Gev are shown in Fig.5. They are 
repulsive and have a range of about 0.3 fm, as one would expect for a short range "nuclear" 
core. The potentials at contact are rather similar in this formalism, and the relative weakness 
of 1=0 scattering is a result of its shorter range. This is an effect of the backward peaking 
diagram D4, which leads to a very short range potential with a large value at contact, and 
carries higher weight in 1=0 scattering. 

Although these Born-equivalent potentials are convenient for use in a meson-baryon 
Schrodinger equation, the actual KN potentials are so strong that they reproduce some 
features of the interaction only qualitatively. For example, the Born diagrams give an 
1=1 scattering length of —0.35 fm, but the Born-equivalent potential (65) for 1=1 in the 
Schrodinger equation for KN leads to a scattering length of only about —0.22 fm. The 
discrepancy is due to higher order effects of Vkn(j") in the Schrodinger equation; we have 
confirmed that the ratio of hfi and Vkn{^) scattering lengths approaches unity in the small- 
a s limit. In a multichannel study one might modify Vkn{ t ) (65) to give the input hfi 
scattering lengths, perhaps through a change in the overall normalization, as a way of 
providing a more realistic potential model of the quark Born amplitudes. 
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V. RESULTS FOR KA, K*N AND K*A; PROSPECTS FOR Z*-MOLECULES 



The channels KA, K*N and K*A are interesting in part because they may support 
molecular bound states if the effective interaction is sufficiently attractive. In contrast the 
low energy KN interaction is repulsive in both isospin states. These "Z*-molecules" would 
appear experimentally as resonances with masses somewhat below the thresholds of m 1.7 
Gev, m 1.85 Gev and ~ 2.1 Gev. Even if there are no bound states, attractive interactions 
will lead to threshold enhancements which might be misidentified as Z* resonances just above 
threshold. 

Plausible binding energies of hadronic molecules can be estimated from the uncertainty 
principle and the minimum separation allowed for distinct hadrons as Eb ~ l/{Mh a d • 
lfm 2 ) ~ 50 Mev. In comparison, the best established molecules or molecule candidates 
have binding energies ranging from 2.2 Mev (the deuteron, which has a repulsive core) 
through 10-30 Mev (the / (975), a (980) and A(1405)). (The / (1710), with a binding 
energy relative to K*K* of about 75 Mev, appears plausible but is a more controversial 
candidate |L3| .) Finally, the state / 2 (1520) seen by the Asterix Crystal Barrel f47 



and Obelix [HS| collaborations in PP annihilation is an obvious candidate for a nonstrange 
vector-vector molecule, with a (poorly determined) binding energy relative to pp threshold 
of perhaps 20 Mev. 

Several candidate Z* resonances which might be meson-baryon molecule states have been 
reported in KN partial wave analyses. The 1986 Particle Data Group compilation [^] (the 
most recent to review the subject of Z* resonances) cited two 1=0 candidates, [Z (1780), | ] 
and [Z (1865), §~] and four 1=1 possibilities, [^(1725), | + ]; [^(1900), § + ]; ^(2150) and 
Zi(2500). However the evidence for these states is not strong, and the PDG argue that the 
standards of proof must be strict in this exotic channel. For this reason these states were 
only given a one star "Evidence weak; could disappear." status. The 1986 PDG also noted 
that "The general prejudice against baryons not made of three quarks and the lack of any 
experimental activity in this area make it likely that it will be another 15 years before the 
issue is decided.". The 1992 PDG compilation []50 makes a similar statement, with "15 
years" revised to "20 years" . 

In their recent analysis of the data Hyslop et al. [||] summarize some previous claims and 
report evidence for "resonancelike structures" [Z (1831),| ]; [Z (1788),| ]; [Zi(1811),| ] 
and [Zi(2074),| ]. The negative parity candidates Z (1788) and Zi(2074) have quantum 
numbers and masses consistent with S-wave K*N and K*A molecules respectively. We would 
not normally expect P-wave molecules; odd-L is required to couple to positive parity KN 
channels, and the centrifical barrier suppresses binding due to these short range forces. 
However, threshold effects which resemble resonances might arise in the full multichannel 
problem, and the very strong spin-orbit force evident in the Poi and P03 KN partial waves 
may be sufficient to induce binding in some channels. A clarification of the status of Z* 
candidates through the determination of experimental amplitudes for the processes KN — > 
K*N, KN — > KA and KN — > K* A in addition to the elastic KN reaction will be an important 
goal of future studies at kaon factories. 

All the S-wave (I,J P ) quantum numbers, in which molecule bound states are a priori 
most likely, £1X6 cLS follows; 
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KA(« 1.6 - 1.7 Gev) : (2, | ) ; (1, \ ) ; 



K*N(» 1.75 - 1.85 Gev) : (1, | ) ; (1, \ ); 



(0, §~) ; (0, i~) ; 



K*A(^2.0-2.1Gev): (2, | ) ; (2, | ) ; (2, ± ) ; 



(1, |") ; (1, |") ; (1, |") 



We can use our detailed model of meson-baryon scattering in the (qs)(qqq) system (q = 
u, d) to identify channels which experience attractive interactions as a result of the color 
hyperfine term. These we again show as weight factors which multiply each of the four 
diagrams D\ . . . _D 4 . Since the overlap integrals these weights multiply are all positive and 
of comparable magnitude, the summed weight can be used as an estimate of the sign and 
relative strength of the interaction in each channel. Positive weights correspond to a repulsive 
interaction. Our results for the hfi "diagram weights" for all KA, K*N and K*A channels in 
(I,S tot ) notation are given below. We also give the numerical values we find for the scattering 
length in each channel given our reference parameter set and masses M K * = 0.895 Gev and 
M A = 1.210 Gev. 



KA(2,-) = - 
v '2 ; 6 



3, -1, +3, -1 



(66) 



a = -0.38 fm . 



KA(l,§) = -iKA(2,§) = l 



3, +1, -3, +1 



a = +0.13 fm 



K*N(1, 



27 



+ 7, +1, -5, 



a = -0.08 fm . 



K*N(1,-) = — 
v '2 ; 54 



+ 26, +5, +2, -3 



(67) 
(68) 
(69) 
(70) 
(71) 
(72) 



a = -0.39 fm 



(73) 



K*N(0,^) = i 



+ 1, +1, +1, 



(74) 
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a = -0.22 fm 



K-N(0,i) = i 



-4, +5, -4, -3 



a = +0.15 fm . 



K*A(2, 



6 -) = - 

2 J 3 



+ 1, -1, -1, +1 



a = +0.14 fm 



K'A(2,?) = 1 



+ 11, -1, -1, -9 



a = -0.20 fm . 



K'A(2,i) = i 



+ 7, +1, +1, +3 



a = -0.86 fm . 



K*A(l,S tot ) = --K*A(2,S tot ) VS tot 



K*A(1,-) = - 
v '2 ; 9 



1, +1, +1, -1 



a = -0.05 fm 



K*A(1,-) = — 
V '2 ; 54 



11, +1, +1, +9 



a = +0.07 fm . 



1 1 

K*A(1, -) = — 
v '2 ; 27 



7, -1, -1, -3 



a = +0.29 fm 
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Evidently attractive forces arise from the OGE spin-spin interaction in the minimum-spin, 
minimum-isospin channels, 



KA: (1,1); 



K*N : (0 



2' 



I). 



K*A: (1,1). 

The two exceptions to this rule are the K*A channels 

K*A: (2,^) 

and 



K*A: (1,|); 

although their weights sum to zero, variations in the detailed overlap integrals lead to at- 
tractive OGE-hyperfme forces in these two channels as well. 

For our reference parameter set we find no molecular bound states; the attractive forces 
are too weak to induce binding. The experimental situation at present is rather confused; 
some references claim evidence for resonances in several channels (see for example 0,^,0 ) , 
whereas other references such as [[30] and |31| conclude that the same phase shifts are non- 



resonant. Our results do not support the most recent claims of resonances [j|], since the 
S-wave quantum numbers of our attractive channels do not correspond to those of the neg- 
ative parity candidates [Z (1788), § ] and [Zi(2074), | ]. However our negative result may 
be an artifact of our approximations, including the neglect of spin-orbit effects and couplings 
between channels. The spin-orbit effects are known from experiment to be very important, 
and might be sufficient to lead to Z*-molecule bound states or strong threshold enhance- 
ments in the attractive channels. Our negative result is based on strong assumptions on the 
form of the interaction; this should be relaxed in future theoretical work, and should not be 
used to argue against experimental searches for possible Z* meson-baryon molecules. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have applied the quark Born diagram formalism to KN scattering. In this 
approach one calculates hadron-hadron scattering amplitudes in the nonrelativistic quark 
potential model assuming that the amplitude is the coherent sum of all OGE interactions 
followed by all allowed quark line exchanges; this is expected to be a useful description of 
reactions which are free of qq annihilation. The model has few parameters, here a s /m 2 q , 
p = m q /m s and the hadron wavefunction parameters, and with Gaussian wavefunctions the 
scattering amplitudes can be derived analytically. The model was previously applied to 1=2 
nix and 1=3/2 Kir scattering with good results. 
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KN scattering is an important test of this approach because it is also annihilation-free 
(at the valence quark level) and the meson and baryon wavefunction parameters and the 
interaction strength are already reasonably well established. Thus there is little freedom 
to adjust parameters. We find good agreement with the experimental low energy 1=0 and 
1=1 phase shifts given standard quark model parameters. (The experimental 1=0 scattering 
length is usually claimed to be very small; we disagree with this interpretation of the data 
and argue in support of a larger value.) A resolution of the disagreements between different 
1=0 KN phase shift analyses, especially at very low energies, is an important task for future 
experimental work. Hyslop |J7| also suggests additional experimental work on the 1=0 
KN system. At higher energies we find that the single Gaussian S-wave phase shifts fall 
with energy more quickly than experiment given standard quark model parameters; we 
attribute most of this effect to departures of the hadron wavefunction from single Gaussians 
at short distances, perhaps in response to the attractive color hyperfine interaction. We have 
confirmed that a smaller hadronic length scale (about 0.7 times the usual nonrelativistic 
quark potential model scale) gives S-wave phase shifts which are in good agreement with 
experiment at all energies. 

We have investigated the possibility of Z*-molecule meson-baryon bound states by ex- 
tending our calculations to all channels allowed for KA, K*N and K*A. Although we do 
find attractive interactions in certain channels, in no case is the corresponding interhadron 
potential sufficiently strong to form a bound state. Of course this result may be an artifact of 
our approximations, in particular the assumption of keeping only the spin-spin color hyper- 
fine term and the single channel approximation. The effect of relaxing these approximations 
would be a very interesting topic for future study. 

There are additional effects in the L>0 KN system which are known to be important ex- 
perimentally, which are not incorporated in our calculations of single channel color hyperfine 
matrix elements. The most important of these is a very large spin-orbit force, which it has 
not been possible to explain as an OGE interaction ]2"2]J2"7|] . Both this spin-orbit interaction 
and the Z* candidates may be strongly affected by coupled channel effects, which we plan to 
investigate in future work. Since much is already known experimentally about the reactions 
KN— >K*N and KN^KA, it should be possible to test predictions of the quark Born dia- 
grams for these channel couplings using existing data sets. Although one might expect OPE 
forces to be important in coupling KN to inelastic channels, such as in 1=1 KN— >K*N, the 
OPE contribution to this process has been found experimentally to be small near threshold 
||. Thus experiment suggests that interquark forces such as OGE and the confining inter- 
action may be more important than meson exchange in coupling KN to inelastic channels. 
We plan to evaluate these offdiagonal couplings in detail in a future study. 
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VII. APPENDIX 



The procedure we advocate for describing hadronic interactions involves simply calcu- 
lating the Born order scattering amplitude for a given process using the constituent quark 
model. In cases where many channels contribute or one wishes to obtain nonperturbative 
information (such as the possible existence of hadronic bound states) then one must extract 
an effective potential (or effective potential matrix in the case of multichannel problems) 
from the Born amplitude and integrate the appropriate Schrodinger equation. The validity 
of this procedure and its relationship to the resonating group method are the subjects of 
this appendix. 

Several theoretical complications arise when considering scattering of composite particles. 
For example, the Hamiltonian may be partitioned in many different ways corresponding 
to different rearrangement channels. Thus if % represents the initial channel consisting of 
hadrons a and b, and / represents the final channel with hadrons c and d, then we may write 

H = H i + V l = H f + V f , (91) 

where 



Hi = --—V 2 R + H a + H b (92) 



and 



Hf = -—V 2 R + H c + H d . (93) 

— * 

Here R is the appropriate interhadron coordinate and [i a b is the reduced mass of the con- 
stituent masses of hadrons a and b. The Hamiltonians H a and Hf, describe the hadronic 
wavefunctions in the initial channel. Thus we have 

H a <fra,n(Q = ea,n0a,n(£a) (94) 

and similarly for the other hadrons. 

In general the wavefunction must be antisymmetrized appropriately and this means 
summing over the various rearrangement channels with the correct weights. Thus we take 
the incoming wavefunction to be a (^a)0b(^b)V ; o(R-) where i/jq is a plane wave and define the 
antisymmetrization operator as 

p 

where P is a quark exchange operator. The Born series for the process ab — > cd may then 
be written as 

(f\4) = ^£R P (-) P 'Wc^o|^ (96) 
PP , \ / 

Here G c is the Green function in a general channel c. 
The Born order expression is 
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(f\T\i) = E(-) P (^0^o|^|P0a^o)- (97) 
p 

This is the "prior" form of the T-matrix. The "post" form uses the potential in the final 
channel rather than the initial channel. The expressions are equivalent if the hadronic 
wavef unctions are exact and the T-matrix is evaluated on energy shell. Note that these 
conditions must also hold if the effective potential matrix is to be Hermitian. 

Little is known about the convergence properties of the Born series. With no exchange the 



conditions for convergence are probably similar to those for simple potential scattering |5T 
Of course exchange scattering is necessarily present when describing hadronic interactions 
from the quark level. At high energy there is evidence that the lowest few Born terms can be 



accurate ||51|| . However, since the potential is strong enough to cause binding in the initial 
and final states, the series will diverge at low energies. Nevertheless the strategy of extracting 
an effective potential can be useful when the Born approximation is not accurate or even 
when it diverges. We shall return to this point below. Despite the theoretical problems, the 
small nuclear binding energy, the small phase shifts seen in K + N scattering, and the lack of 
quark model state mixing evidenced in most of the meson spectrum all suggest the utility 
of the Born approximation. 

It should be stressed that the Born approximation can be useful even when the effective 
potential is very strong. This will be true if the Born term carries information on the 
dominant physics. Then the Born order scattering amplitude may be Fourier transformed to 
yield an effective potential which contains all of the dominant physics and may be integrated 
exactly. Since hadronic interactions must involve constituents, this may not be carried out 
in general, however it will be accurate if the new physics induced at higher order in the Born 
series (such as polarization effects) does not dominate at low energy. As will be discussed 
below, this appears to be true in many cases. 

We now turn to the relationship of this approach to the resonating group method. In the 
following we shall restrict our attention to the single channel case. The resonating group 
Ansatz is then 

* = MU)M&MB.) (98) 

where if) is an unknown function of the interhadron distance. This Ansatz must be antisym- 
metrized. For later convenience, we choose to separate the identity permutation. Thus the 
wavefunction is 

= A* = ^ (l + £ R p pj <p a M. (99) 

Varying the Schrodinger equation with respect to ip and rewriting the resulting expression 
for ip in Lippmann-Schwinger form yields the following equation. 



V>(R) = Vo(R) + 2/i ab J G (R,R')V D (K')^(K')dK' 
- ER P / G (R,R')#;(fa)$(fb) [V| + A; r 2 cl ] p[MQM^)HR')n a ^ h dR' 
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+ 2fX ab £ (-) P J G (R, ROC^a)^^) K(£a, a, R) J P[0a(^)06(eb)^R / )]^a^brfR (100) 



where 
and 



A; r 2 el = 2fi ab (E -e a - e b ) (101) 



Vd(R) = | 0;(ea)0 6 *(a)Vi(ea,a,R)0a(ea)0 b (a)rf«b, (102) 

and G is the Green function for + k^ el . 

The permutation operator in the third term implies that V^+fc^ el is a Hermitian operator 
and may be safely applied to the left. Thus the third term (the kinetic and energy exchange 
kernels) simplifies to 

" E H P / <f>a<f>bP[<f>aM(R)]dUdt; b = f N(K, R)^(R)dR (103) 

where iV(R, R) is the normalization kernel. Because of nontrivial permutation operators 
this expression is damped as R — > oo and hence it does not contribute to scattering. 

We may now iterate Eq. (??) to see that it corresponds to the full Born series ( |96| ) with 
the sum over intermediate states restricted to the appropriate single channel. In particular 
the R —>■ oo limit of the first term in the series corresponds to the Born order T-matrix of 
Eq. (j97D (for the case of elastic scattering). 

Eq. (??) indicates that setting V* = in a resonating group calculation should yield 
a null phase shift. However if one uses approximate hadronic wavefunctions (as is almost 
always the case in resonating group calculations) then a residual spurious phase shift will 
remain. We note that Bender et al. [24j] employ single Gaussian hadronic wavefunctions so 
that one expects small phase shifts upon setting V$ = 0. This is indeed what they found 
for 1=0 K + N scattering. However they obtained rather large phase shifts for the 1=1 case. 
Since the Hamiltonian is independent of the isospin the Gaussian wavefunctions should have 
been equally effective in both cases and one must conclude that there is likely to be an 
error in the 1=1 calculation. Maltman [|52] has concluded that there are indeed errors in the 



hyperfine matrix elements in this reference. 

Solving the single channel resonating group equation is similar to the process of extracting 
an effective potential from the Born scattering amplitude and integrating it exactly. Both 
methods treat the single channel subspace nonperturbatively and hence are successful when 
the single channel approximation is a good one. Both methods fail if (off-channel) virtual 
particles, polarization, wavefunction distortion, and similar effects dominate the low energy 
behavior of the system. As stated above, this does not seem to happen in practice. The 
similarity of hadron scattering amplitudes obtained from the single channel resonating group 
method and from integrating effective potentials has been demonstrated for several cases in 
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